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With the integral representation of Bose functions, the Bose-Einstein condensation of non- 
interacting bosons in a three-dimensional harmonic trap was studied. The relation between the 
particle number and its phase transition temperature was clarified. Some next-order terms in the 
thermodynamic expansions were obtained. We plotted the chemical potential, the mean energy, and 
the specific heat and found most of these properties obtained by using the integral representation 
were almost identical with those of the series representation of Bose functions. 
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After the historical paper by Bose and Einstein in the 1920s, the observation of a Bose-Einstein condensation 
(BEC) was successfully demonstrated in magneto-optical traps (MOT) of dilute alkali gases about a decade ago [J UJ. 
Although the Bose atoms in the MOT are confined and interacting [!, 0, [||, a confined and non-interacting Bose 
system was studied because it is analytic and gives some clue for the lab system when the number of particles is not 
very large. However, the previous studies have been based on series representations of the Bose functions @, 0, [!,[§]• 

In this paper, we study a very fundamental part of the BEC of confined and non-interacting Bose particles by using 
an integral representation of the Bose functions [T3| • We obtain the relation between the number of particles and the 
shift of T c , the temperature of the phase transition, exactly. The chemical potential and the specific heat around T c 
will be displayed to find the phase transition. In principle, we will follow the method introduced by Haugerud et al. 
during the statistical analysis [1, Q . 

Consider a 3D system of N non-interacting bosons confined by a harmonic potential of angular frequency uj. The 
discrete one-body energy eigenvalues of the harmonic oscillator potential are given by e n = H(n x uj x +n y uj y +n z uj z )+eo, 
where n x ,n y , n z — 0, 1, 2, ... and n — {n x , n y , n z }. The zero-point energy is given by Eq = (l/2)h(uj x + uj y + uj z ). The 
distribution function of bosons in the grand canonical ensemble at temperature T is given by 

N = E e ^-,)_ r d) 

n x ,n y ,n z 

where f3 = 1/kgT, ks is the Boltzmann constant, and [i is the chemical potential. For simplicity, we here ignore the 
anisotropy of the trap and consider a typical trap potential of uj x — uj y = u> z = uj = 10 3 sec~ 1 . Since each energy level 
has a degeneracy d(n) = (n + 1)(« + 2)/2, the total number N in Eq. ([1]) is written as 

N = I + V d{n) 
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where z = is the fugacity. 

There are various ways to convert the summation into integrals, and here we use the Euler-Maclaurin's summation- 
5-H ' integral formula 0, H, SEH- Then, the N in Eq. ([2|) turns into the following integral forms: 
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where a = —fifi + /3huj, e 
Bose functions denned by [l(| 



-a _ p 0(^-Hui) 



z e is the effective fugacity [8|, [9J , and F s (a) is the integral representaion of 



Fs(a) 



T{s)J e<*+*-l dX ' 



(4) 



with s > 0. F s is related to the series representation of Bose functions, g s (z) = Y^n=i z n /n s , by way of the effective 
fugacity as F s (a) = g s (z e ). 

The Riemann-zeta function is a special case of the integral as F s (0) = g s (l) = C( s ), an d there are several more 
useful properties 



dF s (a) 
da 



= -F s _ 1 (a), 



(5) 



where s > 1. Note that dg s (z)/dz = g s ^\(z)jz. When a is small, F s (a) can be expanded by using a power series of 
the a as 10 



n=0 

~ r(i- s )a s - 1 + c(s) + ... 

For the first two values of the integer s, F s (a) is given by 

F 2 (a) = C(2)-a(l-lna) + ... , 
F\(a) = — In a + . . . 

Since <C 1, in Eq. © is expanded around — fifj, with the help of Eq. ([5]) as 
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Substituting Eq. ([5]) into Eq. © , we find that the number of particles in the condensate varies with the temperature 
as [USUI 



iV = N 



k B T 

huj 
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In the transition region, fj, — » 0, then, F s becomes the Riemann-zeta function £(s) 
Defining the transition temperature T c where Nq = in Eq. (fit))) , we find that 
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This is the microcanonical critical temperature for the infinite system. A in Eq. (jlip is the one-dimensional inverse 
particle number defined by 

A(N) = C{2) ^ (13) 

A W)- C ( 3 )2/ 3iV l/3 - JVV3- 

Therefore, the ratio of the shift of T c is given by Q 

S(N) = T °~ Tc = . (14) 

{ 1 T 2(1 + A) 1 ' 

If N < 198, then the S{N) is more than 10%. Similarly, N = 2000 corresponds to about 5% and N = 20000 to about 
2.5%. We plotted the S(N) as a function of N in Fig. 1. Since S(N) -> for large N,Tq will be a useful standard 
temperature for discussion. Therefore, Eq. (JTUJ) is written with the definition of To as [g, 

N (T\ 3 F 3 (-pri 3(T\ 2 F 2 (-0n) 

N \T ) C(3) 2 \T ) C(3) 2 / 3 iW3- ^ °J 

The condensate fraction is plotted in Fig. 2 for increasing N The sudden flattening of the curves takes place around 

The negative shift @ in Eq. (fTTj) is expanded up to the order of 7V~ 2 / 3 as 

Tc-To = __4__ 
T 2(1 + A) 

0.7275 1.059 

= -^TJT + JpN + --- (16) 

The chemical potential for N = 200, AT = 2000, and TV = 20000 is plotted in Fig. 3. It approaches 0asT-> T c . Note 
that there is a cross at T — To. 

The specific heat is obtained from the first derivative of its mean energy. The mean energy of the system is the 
internal energy of non-interacting and harmonically confined bosons: 

d(n)e n 



e /3(s n -n) _l ^ ' 



n=Q 

The zero-point mean energy is small in the condensate state. If we apply the Euler-Maclaurin formula Eq. (117[) . we 
can write U (T) in the similar way with Eq. ^ as 

U _ Ne ^ ±m 3 + 3m 2 + \}m + 3 

= iZ£ H 5 — F 4 (a) + — - — Fi(a) H — F 2 (a) 

Tuj GSM* (/3M 3 2((3Fnu) 2 2[ > 
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We may expand U around a = —fl/i. Then, we have 

B = £ + db^'"*' + pb^<-*> + p^"*' + ■ ■ ■ < 2 »» 

Note that the F\ term does not exist here. Replacing ks/h^ with To in Eq. (fT2|) . we obtain the energy per particle 
as a function of T/Tq: 

U-Uo = (T\ 4 F^-fr) /T\ 3 F 3 (-/?M) /rV J2(-/?M) , , 

Nhw \To) ({Z^N-V* \T J C(3) \T J C(3) 2 / 3 A^/ 3 ' 1 J 

We plotted (U — Uo)/Nhuj of Eq. (|2"T|) in Fig. 4. The phase transition is shown at T c . 
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If we take the partial derivatives with respect to temperature in Eq. ([20]) , the specific heat around T c is obtained 
for two ranges: 
(i) When T < T c , 



(ii) When T > T c 



Cy _ (T\ 3 F 4 (-^) [T\ 2 FstM 
Nk B \T ) C(3) \To) C(3) 2 /W3- 



Nk B 



-3T 

The da/dT is obtained from Eq. (fl~2f as 

df = 
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When /i approaches zero, da/dT in Eq. (f24)) becomes zero, too, because F\ diverges. Therefore, C v in Eq. (|22| and 
in Eq. (|23p match at T c . The basic procedure is the same as shown in Ref. |6|]. 

Cy/Nks is plotted in Fig. 5. The specific heat exhibits a sharp jump like the A transition. Since the divergence 
of F\ is extremely slow, it is very difficult to find the continuity at T c in the numerical calculation, but the sudden 
drop of Cy in the vicinity of T c is clear. As previous research pointed out Q, the A-like phase transition in the 
specific heat is seen for any small number of non-interacting and confined bosons. 

The discontinuity of the specific heat at To was obtained as a function of N. If we expand this equation up to terms 
of 1/iV 1 / 3 , we obtain the following form: 



Nk B 
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Except the T c , the point of sudden drop of the slops, the divergence of Fi{—(3jj) is extremely slow. 
T > T c , the N 1 / 3, dominantcs F\ and the approximation is obtained as 



(25) 



Therefore, at 



C v — Cy 
Nk B 



,C(3) 



C(3) 



C(2) - V N 
9.569 



1/3 



= 6.577- 



ATi/3 



(26) 



If N is order of 10 3 , the second term is about 10% of the first term. 

We rewrote the BEC of confined and non-interacting Bose particles in 3D by using the integral representation of 
Bose functions. We obtained the exact relation between the number of particles and the shift of T c and showed how 
the thermodynamic limit of the condensate fraction is achieved around T c with increasing N . The chemical potential 
around T c was plotted for increasing N . It has a common point at To, regardless of the number of particles. We 
plotted the mean energy and the specific heat around T c . Although we used the integral representation of the Bose 
functions, we see that the results are almost identical with previous results obtained by using a series representation 
of Bose functions. 

The author sends many thanks to S. Gnanapragasam and M. P. Das for useful discussions. 
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FIG. 1: Ratio of the shift of the phase transition temperature as a function of N. 
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FIG. 2: Condensate fraction around T c . The solid line is for N = 200, the long-dashed line is for N = 2000, the short-dashed 
line is for N = 20000, and the dotted line is for the iV — > oo limit. 
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FIG. 3: Chemical potential around T c . The solid line is for iV = 200, the long-dashed line is for iV = 2000, and the short-dashed 
line is for N = 20000. 
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FIG. 4: Mean energy around T c . The unit is /iw. The solid line is for N = 200, the long-dashed line is for = 2000, and the 
short-dashed line is for N = 20000. 
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FIG. 5: Specific heat around T c . From left to right, N = 200, N = 2000, and N = 20000. 



